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PREFACE 

This report presents investigations carried out by 
Rao and Associates for the National Aeronautics and Space 
Administration/ Ames Research Center, Moffett Field, Cali- 
fornia under Contract NAS 2-6401. Together with earlier 
reports NASA CR-114576 and NASA CR-2354, this report de- 
scribes the results of analytical studies on fan noise 
conducted under the contract. The helpful suggestions by 
Mr. B. K. Hodder and Mr. D. H. Hickey, the technical 
managers on this contract, are greatly appreciated. 
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f JMMARY 

The method of computing radiated noise from a ducted rotor 
due to inflow distortion and turbulence are examined. Analytical 
investigations include an appropriate description of sources , the 
cut-off conditions imposed on the modal propagation of the pressure 
waves in the annular duct, and reflections at the upstream end of 
the duct. Far field sound pressure levels at blade passing fre- 
quency due to acoustic radiation from a small scale low speed fan 
are computed. Theoretical predictions are in reasonable agreement 
with experimental measurements. 
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LIST OF SYMBOLS 


B number of blades 

V mean relative velocity met by rotor blades 

N rotor revolutions per second 

P far field pressure 

U mean axial inflow velocity 

R distance of field point from center of duct inlet 
'D dipole distribution in the rotor region 
S monopole distribution in the rotor region 

r radial coordinate 

9 angle measured from a reference meridional plane, positive 
in the direction of rotor rotation 

z axial coordinate along the rotor axis measured positive 

from the mid-chord plane in downstream direction 

b blade thickness 

t time 

u axial acoustic particle velocity 

p acoustic pressure 

u c' u l velocit Y fluctuations along and normal to blade chord 

k wave number of the' radiated sound 

A duct length upstream of rotor 

c blade-chord length 

a Q ambient speed of sound 
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l lift distribution along the chord per unit span per 
unit chord length 

e monopole distribution along blade chord 

<p spectral density function 

ft rotor circular frequency 

X specific modal impedence 

V azimuthal angle 

X blade stagger angle 

Y vorticity distribution along the blade chord 

p Q ambient air density 

u circular frequency of the radiated sound 

v circular frequency of convected sinusoidal gust met by 
the rotor blade 

v R reduced frequency of the convected sinusoidal gust 

C distance along the blade chord from mid chord, positive 
toward trailing edge 

n distance perpendicular to the blade chord 

y 

C mv radiation impedence of the open end of the duct 
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Superscripts 


temporal mean 
* amplitude 

+ propagation towards the rotor 

- propagation towards the open end of the duct 
complex conjugate 

i equivalent incompressible plane 

s due to monopole radiation 

d due to dipole radiation 
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1.0 


INTRODUCTION 


Fern noise sources and acoustic radiation from them have been 
the subject of analytical studies at Rao and Associates/ Inc., 
under Contract No. NAS2-6401 with Ames Research Center of the 
National Aeronautics and Space Administration. The results of 
previous studies reported in Ref. (1) and experimental data reported 
in Ref. (2) have well established that disturbances in the inflow to 
the rotor play a primary role in noise generation. The theoretical 
predictions on fan rotor noise presented in Ref. (1) were based 
on radiation to far field from the acoustic sources in the rotor 
region ignoring any duct effects. The present report includes 
the effect of the annular duct enclosing the rotor.. 

I 

In the experiments carried out in the anechoic chamber at 
Ames Research Center on the small scale fan shown in Figure 1 and 
reported in Ref. (2), the annular duct entrance is 14 in. (35.56 cm) 
upstream of the rotor, whereas on the downstream side, the duc\. is 
extended straight through the chamber wall to the exterior of the 
chamber. To take into account the influence of the duct config- 
uration, we extended our analysis to include modal propagation of 
acoustic waves in a duct extending to infinity downstream of the 
rotor but terminated upstream at a finite distance from the rotor, 
as sketched in Fig. 1. 







In our analysis, first we describe: the typical features of 
the nature of the inflow to the ducted rotor in Section 2.. The 
steady state non-uniformity as well as the turbulent velocity 
fluctuations in the inflow give rise to acoustic sources in the 
rotor region as discussed in Section 3. Dipole sources due to 
fluctuating lift forces on the blades and monopole sources due 
to finite thickness of the blades are determined in this chapter. 
Acoustic propagation from these sources is described in Section 4, 
considering the effect of the annular duct on the circumferential 
and radial modes of the propagating pressure waves. The influence 
of the upstream end of the duct on each mode is examined by 
assuming that the duct is terminated by a flange extending trans- 
versely to infinity. Since acoustic energy is transmitted through 
the duct only by the axial component ox acoustic particle velocity, 
we considered only the axial component of particle velocity in the 
pressure waves. Taking into account the reflections of the duct 
modes from the duct entrance, farfield sound radiation is de- 
rived in this Section. Theoretical relations developed in 
Sections 2 to 4 are used in predicting acoustic radiation from 
the small scale fan rotor and the computational results are dis- 
cussed in Section 5. Comparison with experimental measurements on 
the small scale fan in the anechoic chamber at NASA-Ames Research 
Center indicates some improvements that can be made in the analytic 
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predictions. Conclusions drawn from general features of no* 1 
propagation from ducted fan rotors are given in Section 6. 
Recommendations for future effort in theoretical prediction 
of fan noise are also described in this section. 
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2.0 . DESCRIPTION OF TURBULENT FLOW FIELD IN THE DUCT 

The flow in the duct upstream of rotor contains random as well 
as deterministic fluctuations as discussed in Ref. (1). Our present 
interest being in the interaction of inflow turbulence with rotor 
blades, we postulate a suitable description of turbulent velocity 
fluctuations along and perpendicular to the blade chord. Let u and 

V 

denote these velocity components which can be described in terms 
of their mean square values and spectral densities. 

Since at any instant of time, the instantantancous velocity 

fluctuations u (r,6,z.t) and u, (r,0,z,t) must have a unique value 
c ■** 

when the azimuthal angle 0 is increased by any multiple of 2it, vre 
expand u c and u x in a Fourier series in 0 as 


t ' 

P 



u , (r f 0 9 z f t"j! 




; (r,z,t) 


\(r,0,z,tj k^(r,z,t)f 


.e^e.e" 1 ^ 


( 2 . 1 ) 


”here u and u, are complex and random functions with respect to 
c >g *->g 

time t. Now, under the assumption of homogeneous frozen convected 

turbulence, u _ and u. ^ can be written in terms of their Fourier 
c ,g *-,g 

transforms 
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It is important to note here that Eg. (2.2) is only a symbolic 
representation of the random quantities u c and u x . The only physically 
meaningful results are the mean square values u* and u^ . 

The space timt* correlation function R(£,n»CrT) of velocity 
components u^ and u c for frozen convectcd homogeneous turbulence 
is assumed to be of the separable form as 

R(5,n,C,T) - f^ T • f n • (2.3) 

where £,ri,C denote the separation distances along the axial, 
circumferential, and radial directions respectively and x. denotes the 
separation time. Wo further assume that the correlation function 
in the axial and circumferential directions is of exponential form 
while in the radial direction is unity within the correlation length 
and zero outside, that is the "eddy'' is perfectly -orrelated in 
the radial direction. Since the turbulence is frozen convected with 
velocity U in the axial direction, 

f c#T - exp{|£ - U t|/A z > (2.4) 

whore & z is the correlation longth in the axial direction. Since 



the correlation function must be periodic in the circumferential 
direction, we write 


f n ■ ex P{~ In + g 2irr |/£ q) 
g* - • 

where is the correlation length in the circumferential direction. 
Using Poisson Summation formula (Ref. 3 ) , the summation over g 

is transformed into a complex Fourier series to yield 


f 


n 


I 


* 0 /nr 

1 + (gfc 0 /r) 2 


exp (-ign/r) 


(2.5) 


In order to be consistent with the assumption that the correlation 
function is separable, the ratio (SL^/r) has to be considered constant 
with radius. Noting that n = r9, Eg. (2.5) can be written as 


f 


n 


a 1 + <9*/r) 2 


exp(-ig6) 


( 2 . 6 ' 


The correlation function f^ in the radial direction is assumed as 


f 


C 


1 for 0 ^ c < Ay 
0 for c > 


(2.7) 
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Since the correlation function must be unity when t,n,C an a r 
all are zero, the Bq. (2.3) is normalized by dividing it by a 

factor 



U 0 /itr) 

1 + (g* e /r) a 


coth (tr/A 0 ) 


and we obtain 


» exp{|C - U T I A z > 

V ^ e / irr ) 

X tan h {irr/Jl 6 > *Z. — 2 ■ • exp(-ig6) 

1 Qt g l ♦ (gt e /r> # 

( 2 . 8 ) 

Measurement, at inlet of the small scale fan indicated longi- 

t 

tudinal correlation length *> 2 to be same for both components 
u A and u Q . In the absence of specific measurements to the 
contrary we can consider that is also the same, leading us to 
to a reasonable assumption that above Eq. (2.8) gives the corre- 
lation function for both the velocity components u A and u c . 

Now the power spectral density functions (w ) an< ^ ^u A 

c 

can ba evaluated by taking the Fourier transform of the correlation 
function r and we obtain 
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0 1 


iu* I 

c 


c 1 


> 

1 / 



2 I 



lu A / 


*z/irU 

1 + (H w/U \ 2 
z 


X tanh 

g=»oo 


{ur/l 0 ) 


V* r 

1 + (gig /r) 1 


( 2 . 9 ) 


where u£, u* respectively are mean square values of the fluctuations 
along and normal to the blade chord. 


The spectral density of u given above differs from that 

G 

employed in our earlier work reported in Ref. (1), since in our 
present approach, we have assumed the same correlation function 
for u c and u^ components. We also note that, due to the exponential 
form of f^ assumed here, the g th component of above equation is 

different from that employed in Ref. 1. 

1 
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3.0 


ACOUSTIC SOURCES IN THE ROTOR REGION 


The interaction of velocity fluctuations described by 
Eq. (2.2) with rotor blades gives rise to mass fluctuations 
(due to blade thickness) and fluctuating lift forces in the 
rotor region. These mass fluctuations and force fluctuations 
give rise to the so-called "aerodynamically generated noise" 
and can be replaced by corresponding monopole and dipole 
distributions respectively in the same manner as described in 
Ref. (1). 

The velocity field described by Eq. (2.2) is in the duct 
coordinate system, and the velocity fluctuations as met by the 
rotating blade are described in the following section. In 
section 3.2, the influence of these velocities on the blades 
and the resulting distribution of acoustic sources in the rotor 
region are evaluated. The analyses reported here differs from that 
reported in Ref. 1 in that we have included the effect of compress- 
ibility and also considered chord-wise distribution of lift loading 
on the rotor blades. 

3.1 Flow fluctuations met by the rotor blades : 

For the sake of convenience, it assumed that the mid-chords 
of the rotor blades lie in the plane z = 0. Let us denote the blade, 
whose mid-chord occupies the position 0 = 0 at time t » 0, as the 
first blade. Measuring 0 positive in the direction of rotor rotation 
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and counting the blades in the same sense , the mid-chord of the 
jth blade occupies the position 

0 . = 2it^~ at time t = 0 

a b 


Thus the coordinates of a point on the chord of the jth blade 
at a distance £ rearward from its mid-chord at time t can be 
written as 


( V 9 J,C 


* ftt + 2rr 


B 


£ sin 1 


£ cos A) 


(3.1) 


where 


R = circular frequency of the rotor, and 
A *» blade stagger angle 


The leading and trailing edges of the blade section are given 
by £ * -c/2 and c/2 respectively. The coordinates r_ and z 
in Eq. (3.1) are not subscripted by j since they do not depend 
upon the specific blade. 


Since the circumferential velocity flr of v .e blade element 
is independent of time, the fluctuations met by the element 
at £ on the jth blade can be written as 


10 



r 


! 


^ (r C' e o,?' Z C' t, | 


U x (r ,0 ,z ,t)j 

C J,C C 


V- f 1 u „ ,(r f « ) . j(S 

/ I ' J *e ^ g ♦exp{-iu (t- £ cosX/U )} 

J ( u (r,w ) 
a '*»* * t g 1 


X exp{-ig(ftt+2Tr^ - L£i!li) }du> 


In terms of the relative velocity V met by the blade, we can 


write 


C sinl 
r 

C ccsA 
U 


Substituting relations (3.3) in -(3.2) , we obtain the g th 
component of velocity fluctuations as 


, u c,»i _ r r- i c.,‘ r '“ > 


e~^g -expi-i^Tr'i^) 


X exp{-i (g *)+w )(t-y)}du 


I 


r 



3.2 Evaluation of source distributions: 

In determining the effect of fluctuating flow on the rotor 
blades and the resulting monopole and dipole distributions , 
we make the following assumptions: 

(1) Each blade section is considered as an element of an 
isolated airfoil immersed in a "two-dimensional flow" correspond- 
ing to the conditions occur ing at the blade section. 

(2) Only the velocity fluctuations parallel to the chord 
are significant in evaluating monopole distribution due to mass 
flow fluctuations. Similarly, the velocity fluctuations per- 
pendicular to the chord are important in evaluating the dipole 
distribution due to blade lift. 

(3) The flowfield is considered irrotational and inviscid, 
in spite of the fact that the velocity fluctuations u amd u x 
described by Eg. (3.4) do not satisfy the irrotational ity 
condition. It is assumed that the small vorticity in the flow 
does not affect the source distribution for a thin airfoil. (Ref . 4 ). 
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With the above assumptions, the governing equation 
for the perturbation velocity potential due to relative 
'velocity V past the blade is 


\ * + 2“ . + M’ 

a* tt Ct 

o 0 



$ . + <J> 

CC nn 


(3.5) 


where a 0 is the speed of sound; M * V/a Q ; 5 and h are the 
coordinates in direction parallel and perpendicular to the 
blade-chord respectively. 


Applying Lorentz transformation 


£ a ; n « n, t = { t - M?/a Q }/3 ; 0 2 a 1 -M 2 


above equation can be, transformed into the wave 


equation 



• ♦ = 

t t l Z n n 


where <t> « 3 ^ 


(3.6) 


In order to obtain simple closed form expressions for 
the source distributions (monopole and dipole) we follow 
Osborne's approach (Ref . 5 ) which converts the compressible 
flow problem into an equivalent incompressible plane in the 
region close to the airfoil (inner region) using the method 
of matched asymptotic expansions under the assumption 

{e/3 4 } 2 «1 (3-7) 


where p - 


vc 


4 ii n 


o 


I \ 


, v «" <|U+m 
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1 
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This is accomplished by defining a new set of variables in 
the incompressible plane (in the inner region) as follows: 


c/2 3 


uL . <|> i 

2tt3 


<t> 

{ vc 2 /8tt3 3 } 


the equation (3.3 then becomes 


** + 

♦ • t • 

f, Z l l 

? \ n n 


U/e 2 )** 


tV 


which is approximated as 


*\t.t + *\i n i ■ 0 
( ( n n 


(3 JB ) 


under the assumption Eq. ( 3.7 ). 

It should be noted here that the boundary conditions 
applicable for the solution of Eq.(3.8) are those in the 
inner region only. Thus, Eq.(3.8) is required to be solved 
with boundary conditions on the airfoil only. 


The expressions given in Eq. (3.4) for the g th component of 
chordwise and vertical velocity fluctuations, using the above 
transf 9 rmations , become 

i u - - (u >1 . 

e“ i * l 0 •exp{-i.02»ijp ) 


K,* 

i ( 

L* 

( x ,g 

ir | 

J 

' -"00 
X exp{ 


(gfl+ui ) 


u (w ) 

i'9 


X expf-ivNt* - i^/V^Jdw 


(3. 9) 
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where 


v = 2ir; V = 4it|3 2 V/vc 




(4ttBV0cV; 


U c,<T ) 
U i>g r / 


c,0 


f A 

$ * Uc 4 


and v ■ gft+w 


(a) Mass Flow fluctuations and rasulting Monopole Sources: 

A rotor blade subjected to chord-wise velocity fluctuations 
can be represented by a fluctuating distribution of volume 
sources along the blade chord. The volume sources on each blade, 
representing the blade profile, are obtained by solving the 
equivalent incompressible flow problem and thtn transforming 
the result to the compressible case. We therefore first solve 
Eq. (3.8) subject. to the following continuity equation. 


| b(V + u c ^ + 3<t)/3C)“ WSn j 


at n « 0 


(3.11) 


whoro b(r,) is tho thickness distribut ion along tho chord . Tho notation 
f } implies that tho boundary condition is to bo transformed 
into tho incompressible plane. 
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Since we are interested in the noise generated due to 
unsteady fluctuations only, we omit the contribution due to 
steady state velocity V. Also assuming that 

de 

5 -jr « u c the continuity equation (3.11) can be simplified 

to 


Jftl 


n » 0 


1 d Ybfti V „ /0 . r i . i 

jr on -c /2 £ C ^ c /2 


or JA; 

an* 


n t,a 0 


20 ' 


~r ( b ' u c ^ } on ‘° i/2 1 ^ ^ c */ 2 


(3.12) 


From thin airfoil theory for incompressible flow, making 
quasi-steady assumption, we can readily describe the volume 
source distribution along the blade chord as 


% 


8 




lid 

0 * 1 ^ 


(b*u ) 
C,gr 


Since 

i 

8 


4tt0 8 

vc 


eU) 


0 . 13 ) 


tho corresponding nourco distribution in the comprooaiblo 
piano can now bo writ ton an 



i 


u> 


Sc 


(b*u ,) 

c,g 


• (3.14) 


and we note that the source distribution under the assumption 
u Cf g < < a Q » the same as that for an incompressible flow. 
Substituting from Eq. (3.4) for u , we obtain the source distribution 

C 9 V 


g,j, t 


at C on the jth blade as 


9 


.i.%. " L ’ •«" 1 > 

x exp{-i (p-ft+w) (t - “>}•{ + ib^ft/V } dm 


Fox a blade olomont of span dr, wc can consider tho sourco 

s . • dedr concentrated at Q 0. and represent it fc>y a 

3 * J * Q J * 4 

complex Fourier series in 6 . From the source distribution over 

an elemental volume r d6 *,dr • dc located at (r .0 ,C ), we 

o o O 0 0 0 0 

obtain the source strength at r Q due to thickness of the jth 


blade element as 


s g,j lr o )m f u o, j (a ! * • exp{-i ? 2it^i) 

«co 


j - ll 


oc. 


+ “V 2, 
mi' 

X exp{-i(gft + «)(ti*C 0 /V)' } 

X { db/dC Q + i.hgU/V } *dC o dr o d0 o dm 


(3.15) 


In summing the right hand side of Eq. (3. 15) over j from 
1 to U to take into account all the rotor blades , wo note that 






17 



1 



e’xp{*i({7“m')2ir*£^T , “ B 


for {g-"~x r ) « nB with n - - . . - . 0,1 and zero otherwise. 

Consequently, the effect of thickness of the rotor blades on 
noise generation can be examined by considering the influence 
of the monopoles 


$ 


<- 0 > 


B 

2 ir 


00 

rjv 


9 


(w ).»e"’ i fg' 


I 


m 


im '0 
e o 


x oxp { “i (nB H+m ) tt-r.^/V) } 


X { db/dC Q + ibgrfl/V) d ^ o dr o d0 o dw 


(3.16) 
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(b) Fluctuating Blade Loading and Resulting Dipole Sources: 

Due to the fluctuations of the velocity component normal to 
the chord, there is fluctuating lift distribution along the 
blade chord. As in the preceding subsection, we first solve 
the problem in the equivalent incompressible plane and then 
obtain the result in the compressible case, following the analysis 
given by Osborne (Ref. -5) . We introduce a convenient variable 

5 = cos 1 (2r,/c) 

to obtain the following expression for the vorticity distribution 
as a function 6 on the jth blade due to its interactioh with the 
velocity component u ^g* 


Y~ .(6)*sin6 = 

yjJ 


C 


i j U A^ (w ) exp{-i^2ir~i} 


J-l. 


I 


i ! 


# 


■S' V 


* 


X exp (-ivt) « exp (-iv R cos6 M 2 /0 2 ) 

x | [ J o ( V 02) + iJ 1 <V 02) ] 


J o ( V 02) + iJ x <V 02)1 

K ( -iv R /0 2 ) + K (-iv_/0 2 ). C os6 


K ! ^ iV R / 02) + V^V 0 ^ 


V J .m-1 , . , iyn 

/ 1 (cos mfl+ —A sinfl sinmfi) 
m|r 


«n j 




whoro v *» ; and v R » vc/2V 

19 


(3.17) 
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The right hand side of the above equation can also be represented as 


ao 

■/ 


A — 1 7—1 

(w ) e J ‘ y gf»exp{-igr2ir j ~— } *y(6 ,w ) 
X exp { -i (gft+u> ) t } du> 


(3.18) 


where y(6,w ) is a complex function dependent only upon the 
reduced frequency and Mach number . . 

The* lift fluctuation d i ; at (r,0 , ,c> acting on an elemen- 

9 3 0 0 9 s 

tal area dc*dr of the j'th blade can be expressed as 




(3.19) 


The force exerted by the blade on the fluid is equal and 
opposite to that given above , that is 


if gJ “ - e o%.j (?) - dtdr 


(3.20) 


Using complex Fourier series in 0 , the concentrated force 
of Eq. (3.20) can be expressed as 

&f a j P’ V^expt-im'ie-e . r )}«Y *(c) # dcdr (3.21) 

m f 


o 


To include the effect of all the blades, the expression (3.21) 

is summed over J from 1 to 6 and we obtain the force acting on the 

medium due to an elemental volume r d0 dr d£ located at 

oooo 

*$ ,e o' c o* aB 
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1 


d f (r ) 
4 g o 


(u) ) •e~ i< *V y (;.,<•> ) 

r y y 


X exp{-i (nBft+w )t} 


X p{-im C 0 Q/V) , s~^ m ®o»dC 0 dr o d6 o da» - (3.22) 


where g - m = nB, n = -..-1,0,1, 


The resulting acoustic dipole D (r> ) at *• and the force of 

u u o 

df (r ) are related by 

8 ® • 

p o'!t <V - & f a 


which gives & = T 1 df ( r ) 

o 1 o Ka v n O 


o M a o 

K o o y 


(3.23) 


The components of the dipole in the axial and circumferential 
direction are given by 


0 

g,z o 


00 

/ i. sinX* 


•d^ (r ,w ) 

— * *ex'p {-ika 6 1) • dw 

to 
0 0 


V« ( V 


00 

■h 


7“T d ) ’ e X P { -ika -t}.do> 

P 0 ka 6 9 o K o 
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4.0 SOUND PROPAGATION IN JHE DUCT AND RADIATION TO FARFIELD 
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In our earlier work reported in Ref. 1, we considered acoustic 
radiation from the ces in the rotor region to farfield, with- 

out any restrictions imposed by the surrounding duct. The analysis 
presented here is based on modal propagation of pressure waves in 
an annular duct with axial subsonic flow and radiation to farfield 
from the conditions occurring at the upstream end of the duct. 

Tyler and Sofrin (Ref. 6 ) suggest that the modes propagating in an 
infinitely long duct can be considered to yield the conditions 
occur ing at the duct entrance. In the analysis presented here, 
we include the effect of finite distance between the rotor and 
the duct inlet and examine the influence of the impedence of the 
duct opening. In evaluating the latter, we assume that the duct 
terminates in an infinite flange and avoid the complicated 
computations employed by Lansing in Ref. (7 ) for the case of 
unflanged duct. 

The modal propagation of acoustic waves due to an isolated 
source in an infinitely long duct is described in section 4.1, 
which is followed by a discussion of the integrated effect of the 
rotor-generated source distribution in 4.2. The effect of 

v. 

duct termination at a finite distance upstream of the rotor 
and acoustic particle velocities occuring at the duct-end are 
discussed in section 4.3. The method of calculating radiation 
to farfield and spectral density of farfield sound pressure level 

are described in section 4.4. 
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4.1 Pressure waves in an infinitely long annular duct due to 
source distribution in the rotor region. 


Let us consider sure waves propagating in an infinitely 
long annular duct due to monopole and dipole distributions gxven 
by Eq. (3*16) and (3.27) respectively. Following the analysis 
given in the Appendix , the pressure at r due to an elemental 
monopole S (r Q ) is given by 


Ap® \r) 


00 

(Bka o 0 o /2*iy f ) e" 1 ^ £ e ' im ' 9 0 exp Uka^/V) 

_ m ' 


xY y — - 

o 4,1 V* 


im(0-0 ) 


Tea 
o 


0 % (r > E e„‘ r o> 


X exptik _ (z-z )> • exp(-ika t) 

Z 0 0 


x 


Sr dc dr d0 dw 
0 0 0 


(4.1) 


OftirnGxr n 0F TTO 

0WGJ2VAL PAG£ is P00ft 
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and that due to dipole (r» o ) is given by 


AE V (r) 


B <>o 


00 

r f g <“ e 

m ' 


-im'0 


X exp(~im'C 0 fi/V) 


*1 I 

m~*- oa y=0 


' m cos 1 

-6 X, 

! o o 


+ k± sinX | 


0 


B im(9-9 0 ) (r)( ) 
my my o 


X exp{ik± (z-z q ) ) * exp (-ika Q t) 

X d ^ 0 dr o d0 o dw (4.2) 

where k± «* — (±a - M ) , ka = (nBft+u ) 

z ft 2 o o , 

p o 

and ct = { 1 — 2 ^ 

In equations (4.1) and (4.2) , upper sign is chosen for down- 
stream propagation and the lower sign is chosen for upstream 
propagation . 


The effect of all the sources in the rotor region is accounted 
for by integrating the right-hand side of Eq. (4.1) and (4.2) 
between the 1 imits 


r h i r o i' r t 


0 < 9 < 2n 

— o — 

-c/2< C 0 < c/2 


where c ■ z /cosX 
0 0 
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In carrying out the integration over 0 Q , we note that 


exp{ -i (m+m ' ) 0 }d0 
0 0 


2ir if m=-m r 


0 otherwise 


Thus the myth mode of the duct pressure at (r,0,z) due to 
monopoles in the rotor region can be written as 


CD 

1 (r) » BP o a o f * 

r ,m u ( 4uaA J J u c . 


my 1 r u 

h 


(w )e“ i ^a*E (r ) l. 
g v my' o' b 


X E(r) • exp{ im0 + ik z - ika t> dr du 


where 


(4.3) 


C/ ^ 

= J” exp{iC Q (ka Q /V - k* cosX))X{g|-+ ibgfl/V} d? ( 


-c/2 

and that due to dipoles as 


(4.4) 


P 0,my 


Bp o f r rt - i4 

J { u a,/“ )Ve s 

B 0 my -oo r h 


X (7^f X+ sin\>E mlJ (r 0 ) h 
B o r o 


x * exp{im0 + ik^z - ika Q t) dr Q d<*) 


where 


(4.5) 


w )»exp{ic 0 (^ - k* cosX)}dC 0 


-c/2 
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The Equations (4. 3) and (4.5) can be written as 


•my 


00 

/ 


A (w ) exp{i (m0-ka t + k'z) }E (r)du> 
my 17 o z my 


(4.7) 


where 

A (w ) * 
m\s 


B p« a « 
0 0 


4naA 


my 


t 

/ u e (r )I.dr 

c,g a my o b o 


(4.7a) 


for acoustic propagation due to monopole distribution in the 
rotor region, and 


A w <“ > 


BP. 


A 


a kct . 

4ir -g- A 

o n 


sinx) 

oo 


(4.7b) 


X E (r )I„dr 
A my o l o 


i / 


for acoustic propagation due to dipole distribution in the 
rotor region. The superscripts + and “ respectively indicate 
wave propagation in the downstream and upstream regions of the 
sources . 

The speed with which a mode propagates in the axial direction 
is given by a Q k/k“. Since k z depends upon a and M q as given in 
Kq. (4.2) , let us consider the various conditions governing 
wave propagation or decay. 

a) mode wave number k < k/8 : 

m \ j o 

The value of a for this case is real and positive, leading 
to real values for k*. Consequently, Eq. (4.7 ) indicates 
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that the inodes propagate upstream without decay 
for all real and positive values of a. However, 
the propagation in the downstream direction takes 
place only when a > M q . 

b) mode wave number k = lc/3 Q s 

The value of a. for this case is zero, and we find that 
the mode propagates only in the upstream direction with 
a speed of (a/M Q ) (1-M 2 ) . 

c) mode wave number k > k/0 : 

m\t o 

For this case we can write 

k^ «=» — i« - M ) . 

V 

with a = |l-(0 o k my /k) 2 | ,s 

We find that the mode propagates only in the upstream 
direction. However, the amplitudes of pressure and 
axial perturbation velocity in the mode would decay as 
exp {kaz/0 2 } , where z is the distance from the rotor 
and is negative in the upstream region. Since a 
increases with the mode order , we note that higher 
order modes decay at a faster rate. 
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4. 2 Effect of duct termination upstream of rotor. 

Let us now consider the situation that the duct has infinite 
length on the downstream side but is terminated at a location 
z - upstream. The myth mode of pressure wave in the upstream 
direction due to sources in the rotor region is given by 




(to )«exp{i(m0 - ka t + k”z)}»E (r)dto 

0 Z ffl\X 


oo 

J 


p (to )-dto 

*my 


(4.8) 


where A ffl ^ (to ) corresponds to superscript (-) in Eqs. (4.7a) and 
(4.7b ). Depending upon the impedence of the open end, the pressure 
waves will be reflected and travel back giving rise to right- 
running waves in this region upstream of the rotor. The pressure 
in each myth mode of such reflected waves can be described by 


f 


p my 



A* (w ) »exp{iin0-ika t + ik* 
my r o z 


(z+A)}^ (r)du> 

fflji 




(to )dto 


(4.9) 


i 

where A is a complex constant which depends upon A , l, and the 
impedence of the open end. The superscript (+) denotes that the 
above equation applies to the reflected right-running waves. We 
note that z is measured from the rotor, whereas the origin of these 
right-running waves is at the open end given by z ■ -A, 
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The axial acoustic perturbation velocity u~ due to pressure 
waves given by (4.8 ) or (4.9 ) is evaluated from the linearized 
momentum equation in z direction 


& + » 


7z u± 


1 3 ± 
0 o 7zP 


(4.10) 


where O is the mean axial velocity of the stream. Using Eqs. 
(4.8 ) and (4.9 ) for p* we obtain 




k 4 

z 

O A 


4 1m > 

"o a o 



p ("‘ ) 
1 ml I 


P a 
H o o 


with x 


k-M k 



(4.11) 


The factor x 


in termed an the "specific modal impodonco ratio." 


The perturbation velocities in the circumferential and 
radial directions are not considered since they do not carry 
acoustic energy across a transverse section of the duct. 

To determine the pressure and velocity at the duct termination 
at z ■ -Jt, resulting from an myth mode of left running pressure 
wave, we follow the method suggested by Morse in Ref 8 and elabor- 
ated by Kaji in Ref. 9. We assume that the duct is terminated 
at z ■ -S, by an infinite transverse flange, but extends to infinite 
length in the positive z direction as shown in Fig. 2. Considering 
a single myth mode of fluctuating distribution of axial velocity 
.of the type 


29 



1 


U mu * %/’V r)-e ='P { '- i(ka o t * b6 > ) 


over the annulus r. < r < r , , 

n — — t 


it is shown in Ref. 9 that the resulting pressure distribution 

over the annulus can contain other radial modes tha.i E . (r) . 

my 

In terms of modal representation, we obtain 


P»p- -°o a o u mp -eJt e ( - i(ka o t ‘ ” e)) 
00 

v V E (r) 

A '’mv m\> ' 


(4.12) 


The subscript v indicates all the radial modes resulting from 
the radial mode y describing . The open-end radiation impedence 
ratio is discussed in detail in Ref. 9, and is given by 


r V „ ik 


where A 


and I (t) « 
my 


f I (T) I (T 
L. / mu mv 

IV 0 J (t 2 -k 2 ) ^ 

/* 

J {E mp tr))!-rdr 
h 

J J m (Tr), 1 V r) 


J (tr)-E (r) *rdr 
m my 


(4.13) 


Based on coupling impedences c^computed for a plane wave by Morse 
in Ref. 8, we assume 


<< 

mv my 


to simplify Eg. (4.12) to 

b « “0 a »u »C^ 
*my K o o my '’my 


(4.14) 
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Noting that the pressures and velocities at z - -l in the above 
equation are obtained by summing the left-running and right-running 
waves, we obtain 


P»U ( “ > * 


U 


> > + U m,> > 


P (ti) ) 

= - P a c w 

+ * « o o m\i 


(4.15) 


mp ' ’ ”mp 

Substituting for P^y( wf ) from Eq. (4.11) into the above/ we get 

_P 

(4.16) 


ul (to ) 

TnUi 


± 

^mp 

x- + ? 


u (to ) 

;«p ' 


x + + ci 


HE 


Using Eq. (4.11) we can also obtain 

4 ( “ > x+l w < “ > 

— - - 5 -= (4.17) 

Pm U ( “ > x V" > 

We note that v ± and c 11 are functions of w . 

m\i 

At high frequency i.e., when k ■* »: 

a -*-.1, ' k~ ■+ — (±1-M ) 

Z a 2 0 

e o 

X* *1 and * 1 

then we note from Eq. (4.16) and (4.17) that both 
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u w ( “ > 
u mu ( “ > 


and 


p + (w ) 
^my 

p~ (to ) 
*my ' 


which means that the myth mode passes through the open end 

without reflection. However, at frequencies near cut-off i.e. 

when k ■*> B k , from above: 
o n?y 


o. 


kM 

0 


7 and x± k 


n 


then we note from Eq. (4.16) and (4.17) that both 


u + (to ) p* (to ) 

— and P J!sU. 


U mu ( “ > 


P (o> ) 
*my 


- - 1 


which moann that the myth motion near cut-off do not paint through 
the open ond and nro coipo lately reflected bnck into the duct. 

At I ul hi mod ! tt I n I lie value of plnyw an 

imporl nnl role In defining l ho com 1 1 1, ionn nt thn duct entrance 
on which wo can onl (male the m’ounl lo radiation to farflold. 
l**o r example, tho velocity f luctuutulon at •/, •• -9, in the myth 
mode given by 


u ,,{(0) « u“(<o) + u* <to) 
wy w?y wy 


reduces to 
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using Eqs. 


X* - X- P„ u <<■» 

+ ix " - 

y + pax 

* *my K o o* 
(4.11) and (4„16). 


ubstituting for p (w) fror.. Eg. (4. 8 ), we obtain 


u m u ( “> 


w 

■f 


%<“> X - X 


j pax” x + + C** 

_oo o* A *my 


X exp{i(m6 - ka A t - k~A))*E (r) »du) 

" o z my 


(4 


We note that velocity u in the myth mode travelling upstream 
to the duct entrance together with reflection due to impedence 
t’my' gives rise to the resultant velocity u given in the 
above equation. 


.18) 
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4. 3 Far Field Radiation from the Duct Opening. 

In the preceeding subsection, we calculated the myth mode 
of axial acoustic velocity at the duct opening due to sources in 
the rotor region taking into account the impedence of the duct 
opening. To calculate the far-field pressure at a point ( R , HO 
as shown in Fig. 2, we follow the method suggested by Tyler and 
Sofrin in Ref. 6. The radiation at (/?, HO is considered to be 
that from distributed monopoles of strength 2u r d6 dr. the 

Will 

factor of 2 being due to the fact that radiation only to the outside 
of the flanged duct is considered. Tyler and Sofrin ignored 
reflections at the duct opening and used the velocity in the 
upstream travelling waves given by Eq. (4.11) to estimate these 
monopoles. In the following analysis, we included the duct open-end 
impedence and the resultant velocity given by Eq. (*.18'.' 

The farfield pressure at a point (£,¥) due to a source of 
strength 2u^r d0 dr and curcular frequency ka Q located at 
(-i,0,r) is given by 

r ikh 

» > - / - ik «o°c. l2 % ( “ )r d9dt >-rar du (4.19) 

From Fig. 2, we note that tho distance h from the source is 
given by 

h = {/?*+ r J - 2ri?sinty cos©}* 4 

« 

i 

2 R * r cos0 sinV, for small values oc r/R. ' 
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The pressure at (/? , V) due to all the distributed monoples 
at the duct opening is given by integrating (4.19) over the area 
of the duct and we obtain 


W' 


r 2 tt 

*» - / { I - ik *o0o- e 

too 11 h *D 


ikh 


4trh 


2u ,, (w) r d0 dr dw 


my 


(4.20) 


Substituting the value of h in the above equation, and ignoring 
higher order terms we get 




(», D - J ika 0 P 0 >- 


u „u <"> 

2tr/? 


- r h 0 


X exp{ikR - ikr sin ¥ cos0}«r d0 dr dw 


(4.21) 


Substituting Eq. (4.18) for u^ into the above Equation and 
noting that 


r 


2 IT 


J exp{i(m0 - kr sinY cos0)} d0 ■ (-i) m 2xJ m (kr sin H) 


•6 


the Equation (4.21) becomes 


my 


- / 

•'-c 


(-i) m + 1 -V— • A” (w) »I (k sin *40 

« * + C. p 


X exp( i (ki? - ka t - k” i)> dw 

Q) 2 


where I , (k sinV) 
m\x 


& 


(kr sin'D’E (r) r dr, 
m m\i 


u 

a notation previously employed in evaluating c* . 


* (4.22) 
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08 (v) 


{f F 8 } . {f F 8 } * <“> 

c iy 


due to monopole distribution in the rotor region and 


* p d (v) 




due to dipole distribution in the rotor region. 


Where 


(-i) m+1 k 


+ . u 
Y + C 
* ^mu 


I (k sin Y) , 
my 


Bp a 

> S a 0 £ 

4iroA 


/'• 


E (r ) I. dr 
my o o o 
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RESULTS OF COMPUTATIONS AND COMPARISON WITfl 
EXPERIMENTAL MEASUREMENTS 


Based on modal propagation within the duct and reflections 
at the open end, rad i-'t- ion to far field from the small scale 
fan is computed using the equations derived in the preceding 
sections. The small scale fan configuration is sketched in 
Fig.l' and its detail description is given in Ref. 1. The 
results of such computations are discussed in this section, and 
are compared with measurements. 


In our acoustic computations, the fan data and other physical 
quantities used are as follows: 


§ 


a 

a 

s 

It 

8 


i ) 


B * 15, r K ■ 0.249 ft., 
” (7.59 cm) 


0.598 ft., 
(18.23 cm) 


Blade chord (c) 

Stagger angle (X) 

U «83.5 ft/nec, 
(25.45 m/fioc) 


hub 

0.125 ft. 
(3.81 cm) 

48° 


‘t 

pitch 

0.117 ft. 
(3.57 cm) 

63° 


N «125.3 rps 


tip 

0.115 ft. 
(3.51 cm) 

69.8° 


p « 0.002378 slug/c. ft* 


a ■ 1128.6 f i /roc 
0 


(0.001 226 grnm/c. cm) (344 m/nor:) 

The support spider located at 7 in. (17.78 cm) downstream of 

the rotor has 8 vanes of 1.625 in. (4.13 cm) chord and 0.25 in. 
(6.35 mm) thickness. 

Restricting our attention to noise luvols at blade panning 
frequency, we have 


k - [fi. 

a o 


2nBN 

a. 


10. 6 per ft. (34.33 per motor) 
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From upstream mode cut-off condition 


k > *o k m>i' 

we find that only the following mu modes need be considered i 
m —4 , . . • • —1 |0|l|o«t*«4 

U - 0 

From the relation g+m - B, we note that each gth component gives 

rise to the corresponding mth mode. Consequently, the values of g 

need to be considered are between 11 and 19. For each of the 

(mu) modes listed we computed the radiation impedence of the 

my 

open end of the annular duct, using the method described in the 
preceding section. The real and imaginary parts of 5^ thus 
computed are shown in Fig. 3 

In the experiments carried out on Rao Fan-2 in the anechoic 
chamber at NASA-Ames Research Center, there are no discernable 
upstream disturbances except for turbulent fluctuations in the 
inflow. Consequently, we consider the random moncpole and dipole 
distributions due to blade thickness and blade loading resulting 
from the interaction of the random velocity fluctuations with the 
rotor. 


( i 


The turbulence measurements taken at the fan inlet (private 
communication from Mr* B. K. Hodder) indicated that the mean square 
value of velocity fluctuations in the axial direction was 0.75% 
of the axial velocity of 83.5 ft/sec. (25.45 m/sec) Nearly the 


■ft 
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same intensity was measured for the velocity fluctuations along 
the circumferential direction. Consequently, for our acoustic 
computations we used 

« 

{u^-Cu*}* a 0.6263 ft/sec. (0.1908 m/sec) 

Based on measured auto-correlation function for axial velocity 

fluctuations and a "frozen-convected" assumption with a convection 

velocity of 83.5 ft/sec (25.45 m/sec) , we obtained 9. » 25 in. 

z 

(0.635 m) for us in our computations of the turbulence spectrum. 

A value of « 3 in. (7.62 cm) is assumed to conform with our 

earlier computations presented in Ref. 1. The spectra of the 

turbulent fluctuations, ^u, „ and are then computed from 

i.,g c,g 

Eq. (2,9), for substitution in Eqs. (4.25) and 4*26). The mean 
square value of blade lift distribution computed from Eqs. (3*17) 
and (3.19) for ui*0, is shown in Fig. 4, 

Xn the preceding section, we cttricusncd the left and right 
running sound waves resulting from the acoustic sources in the 
rotor region, which in turn depend upon the random velocity 
fluctuations met by the blades. To examine the effect of the 
particle velocities in the pressure waves on the generation of 
the acoustic sources, we carried out the following calculations. 
Based only on u (i.e., inflow turbulence) we obtained u 
from Eqs. (4.8 ) and (4.11), and then from Eq> (4 . 16) for 

each value of g and corresponding mode (mt). Since thete velocity 
components are random, one can compare only their spectral 
densities, which depend on radius due to the function E my (r) in 
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in Eq. (4.7,). From the behavior of E (r) for each of the myth 

my 

modes considered as shown in Fiq. 5, we note that the velocities 

u* y are largest at tip radius. The spectral density of 0 u 

also is a function of radius as can be seen from Eq. (2.9 ). 

Consequently, we calculated 0 - , 0 U + and 0 at w = 0, 

my my u x,g 

and r • r.. and listed them in Table 1. We observe that 0 

t u. 

A ' g 

dominates and similar situation can be expected, when we consider 

the outcome of monopole sources due to the velocity fluctuations 

u c,g* Consequently, we can ignore the interaction of propagating 

sound waves with the blade row, and justify the use of Eqs. (4.25) 

and (4.26) with the corresponding values of 0 and 0 

c,g u ±.g 

From the blade loading based on the random fluctuations of 
u x , we computed sound pressure spectrum according to Eq. (4.26) 
at distance R « 7 ft (2.13 m) and angle H* ranging from 0 to tt/2. 
To reduce the amount of computational work, we restricted our 
attention to the spectral density at the blade passing frequency 
by setting u « 0 and v a 2 ttbN in Eq. (4.26). The sound pressure 
levels thus computed are shown in Fig. 6. 

In computing noise from blade thickness effect according to 
Eq. (4.25), we used a blade profile described by its thickness 
distribution b(c) given as 

\ b(c) ■ constant • (1 + j-£) {-log(l + §£•) + log 2)i 

G G C 

with “c/2 < ( < j 
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By choosing 0.27182 for the constant in the above equation, we 
obtain a profile with a rounded leading edge and a maximum thick- 
ness of 10% of chord occurring at 37% from leading edge. The 
source-sink distribution e(c) is obtained by substituting 

b(C) from the above equation and u from Eq. (3.4) into : 

c ,g 

Eq. (3.14). The results of our computations for farfield noise 
due to blade thickness effect are also shown in Fig. 6. 

The SPL dB levels shown in Fig. 6 for noise from the blade 

loading and blade thickness are from the spectral density per 

Hz of the acoustic pressure at blade passing frequency. Since 

the fluctuations of u and u are considered random, one can 

x c 

add the mean square values of sound pressures from these two 
contributions to obtain the combined effect. However, we observe 
that the noise contributed from blade thickness effect is neg- 
ligible in comparison to that from blade loading* 
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The experimental data given in Ref. 2 show the spectral analysis 
obtained using 25 Hz band-width. From our previous computati >ns 
of sharp peaks in sound pressure around the blade passing 
frequency, as re? nd in Ref. (1), we can assume that integration 
over a 25 Hz band -v/iuch leads to a spectral density 12.1 dB higher 
than those computed from Eq. (4.25) or Eq. (4.26) at the blade 
passing frequency. The sound pressure levels thus estimated are 
shown in Fig .6 along with the measured data. Apart from the 
noise level, which is proportional to turbulence intensity, 
the theoretical estimates do not show, the strong directivity 
found in the measured data. 

We also computed the acoustic radiation to farfield from the 
sources in the rotor region, completely ignoring the effect of 
the duct, for the sake of comparison. The sound pressure levels 
thus estimated are shown in Fig. 7 and we note that they do shew 
the type of directivity observed in the experimental data. 

The inability of our theory in predicting the directivity 
of noise from a ducted fan can be due to any or all of the 
following reasons: 

(a) The effect of the confining duct and its finite length 

♦ 

are examined in this report by considering the radiation impedence 
of an annular duct terminated by an infinite flange, as sketched 
in Fig. 2. The duct configuration used in the tests, shown in 
Fig. 1, is a cylindrical section upstream of the rotor, and is 
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annular downstream. The propagating modes, and the open-end 
impedence for each of these modes in such a duct configuration 
may be different from those considered in this report. Modelling 
the bell-mouth configuration of the fan inlet by an infinite 
flange can be yet a suitable method. Even though the computations 
for a plane wave (i.e., m =* y = 0 mode) presented on page 472 of 
Ref. 10 show that the absence of the flange does not appr ciably 
alter the radiation impedence of the duct opening, there is no 

indication that similar behavior exists for higher order modes. 

** * 

(b) An ( my )th mode incident on the open end of the duct can 
give rise to various other radial modes. These (m,v)th modes 
irrespective of duct cut off conditions can affect radiation to 
far field. In our analysis we ignored all such cross-coupling 

of modes. Computations for a plane wave presented on page 336 in 
Ref. 8 show that there is negligible cross-coupling at high 
frequencies. However, there is no indication that cross-coupling 
can be assumed negligible for higher order modes. 

(c) The strong directivity occur ing in radiation from blade 

loading directly to farfield (as shown in Pig. 6 ) appears to be 

due to the dipole nature of the sources. In evaluating the farfield 
radiation from conditions occuring at the duct entrance, we con- 
sidered only the axial component of the acoustic particle velocity 
at the duct termination and the corresponding monopole sources. 

We have ignored the gradients of velocity, and thus lost the dipole 
. character of the sources. In particular, we have not included 

the tangontiul components of the particle velocity and its gradient 
in the conditions occuring nl. Hn> ilinM-. opening. 
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1. Consideration of chord-wise distribution of the blade 
lift loading is important in evaluating sound at wave lengths 
comparable to the blade chord. 

2. The farfield noise, predicted under the assumption of an 
annular duct terminated by an infinite flange, does not show 
considerable directivity. On the contrary, experimental 
measurements indicate about 10 dB drop on the sideline from 
that on the axis. 

- 3. Sound pressure levels computed for an axis point including 

the effect of finite duct length upstream of the rotor is 
comparable to those computed for a free rotor in the absence 
1 of a duct. 

I 4. The acoustic particle velocities in pressure waves reflected 

at the duct opening are negligible compared to the inflow 
turbulent velocity fluctuations based upon which pressure waves 
propagating upstream towards the duct opening are calculated, 
consequently, the second order effect on sourc< 3 in the rotor 
region due to interaction of the propagating pressure waves 
! with rotor blades can be ignored. 

i' 
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Comparison of theoretical estimates with experimental measurements 
and a critical examination of the assumptions employed in our 
analysis lead us to the following recommendations for future 
effort in fan noise evaluation and reduction of far field noise. 


1. Space-time correlation measurements of flow in the duct 
are of prime importance in making theoretical predictions of 
noise from the ducted rotor. The spectral density distribution 
of fluctuations of u c and u^ depend upon the measurements of 

and the autocorrelation function. Even small steady-state 
velocity deficit of 0.5 ft/sec (0.152 m/sec) when occurring over 
a small region of the fan inlet, can give rise to on axis-noise 
comparable to that presently estimated from the inflow turbulence. 
Detailed measurements at fan inlet, as suggested here, can detect 
such velocity distortions. 

2. The assumption of an infinite flange at the opening of an 
annular duct is employed in our analysis to avoid complexity 

j and excessive computational effort. Also, the "cross-coupling" 

of radial modes at the duct opening is assumed negligible. The 
| discrepancy between the predicted and experimentally determined 

j directivity patterns may be due to the above simplifying 

,) 

assumptions, and also partly due to some other noise sources 
overlooked in our analysis. Even though our primary interest 
j is in reduction at its sources, modal transmission of acoustic 
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waves in ducts with cross-sectional discontinuity similar to 
that occurring in practice , and reflection of higher order 
modes at the duct opening requires further examination to 
formulate reliaole fan noise prediction methods. 
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APPENDIX 


Pressure waves in an infinitely long annular duct due to 

an isolated source. 


Let' us consider a. mono pole source 
S (t) - £'*e' rikfc o t 

f * 


(Al) 


inside an annular duct r^ < r < r fc with an axial inf low velocity U a « 


The governing oquntion for pressure is the acoustic 
equation which has undergone a Galilean transformation since 
we have a moving fluid of subsonic Mach number M Q = U ./a Q . 


1 


v 2 p + 2“oi!E + M 2 i!£ . 
a Q 3z3t 0 3z 2 



6(r-r 0 ) 6(6-0 o ) 6 (z-z q ) 


where V* ■ 


9 2 + 1 3 _ + 1 _ 9 _^__ + 9 ^_ 

3r 2 r dr r 2 3r 2 3^ 2 


(A2) 


The right hand side of Equation (a 2 ) represents the 
Lighthill source term. The boundary conditions on pressure are 

If - 0 at r - r h and r t 
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Using the Lorentz transformation 

r ® R; 0 “0/ z - BZ, t » (T -M Z/a )/B 

0 0 0 

where 0 2 - 1-M 2 , the Eq(4.2) can be written as 

o 0 

1_ 3 2 P _ (a 2 . 3 t 1 3 2 + 3- 2 - ]p 

a 2 3T 2 (3R 2 3R R 2 39 2 3Z 2 j 

- -i£a Q P 0 5*|5 (R-R o ) 6(0-O o ) fi(Z-Z Q ) • exp{-ika o (T-M 0 Z/a 0 )/^ 

(A3) 

and the boundary conditions become 
- 0 at R - r h and r fc . 

i 

We sook tho solution of Eq. (A3) in tho form 
k 

p s " -ig a 0 P 0 s g(<V^ 0 > • exp{ -ika Q T/B } (A4) 

where the Green's function g (/?//? Q ) is defined in the transformed 
coordinates . 


Substituting Eq. (A4) into Eq. (A3), we find that g { R/F Q ) 
satisfies the differential equation 

3 2 , 1 9 , ia 2 | 3 2 kl 

3R 2 R 3R R 2 3 0 2 3Z 2 0 2 


9 (R/R 0 i 


- -6 (R-R ) 6 (0-e ) 6 (Z-Z ) * expfikM Z/0 } 
0 0 0 o u 

and the boundary conditions 
- 0 at r . r h and r t . 


(A5) 


(A6) 
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1 


? 

t 


1 


T 


Solving Eq. (A5) subject to the boundary conditions ( A6 ) 
by the method of separation of variables, we obtain 

oo oo 

« <V V - I I 

4TT - x A 

m=-» y=o B 0 m\i 


X ^ui (R)1 W R o >c ex P (ikM o Z o /0 o ) 

X exp{±i-^ (z-z o )l ( a7) 

where 


In Eq. (a 7) we take + sign for the region downstream (on 
right) of source location while for the region upstream of source 
location (on left),, the negative sign is chosen. 


Rewriting in original variables, the Eq. ( A7 ) becomes 


<r/ v -i i tt 


e im(0 - 9 6> 


lrrt 

m»-« 0 y=0 4ir-r- A 

P MM 
0 


i ) 


x E my (r ] ): E my (r o ) * ex P (ikM o V * exp{ 11 ka ( z *’ z o * / K } 
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(A0) 


and the acoustic pressure at (r,8,z) due to a monopole is given by 


p s (r,e,z) - "ijfVo i'srfr/O 


X exi t-ika Q t) • exp (-ikM o z/0 2 ) (A9) 

Let us now consider the acoustic pressure at (r,6,z) due to 

A 

a dipole £ (t) = ^e’ ika o t at (r,0,z), if being the strength of 
the vector dipole and ka Q * (nBft+w * ) . 

The acoustic pressure due to dipole is easily determined by 
differentiating Eg. (4.4) in the direction of dipole axis and re- 
placing the monopole strength by dipole strength. Thus we obtain 


p d (R) = 'ijra 0 P 0 D * v 0 g(«A 0 )exp. (-ika 0 T/8 0 ) 
0 


(A10) 


where 


A ft A a A a 

7 o * 1 5E + kjiT 0 +k 5^ 


Writing in original variables, Ej. (4.10) becomes 


p d (r)-<-i£ a p { D . | SL + d • - + D • $ |S— } 

p 1 1 ]L<o M o 1 r 3r 6 r~Je~ z p d <Fz7 


o o 


X exp(-ika t) x exp(-ikM z/3 2 ) 


(All) 


AAA 


where (0 . D a , D ) are the components of $ in r,0 and z 
r o z 

directions respectively. 
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Table 


Mode 

my 


«(- 4 , 0 ) 

(- 3 , 0 ) 

(- 2 , 0 ) 

(- 1 , 0 ) 

( 0 , 0 ) 

( 1 , 0 ) 

( 2 , 0 ) 

( 3 . 0 ) 

( 4 . 0 ) 


Z. Spectral intensities of inflow turbulent velocities 
and the particle acoustic velocities in the 
propagating duct-nodes* 


0 

u / 

-i,(7 

0 _ 
u 

my 


$U + 

' my 

(turbulence) 

(incident on 

(reflected from 



duct opening) 

duct 

opening ) 

.430 

(10" 5 ) 

1.121 

(lO” 8 ) 

.068 

(lO” 8 ) 

.480 

(10" 5 ) 

1.464 

(lO” 8 ) 

.227 

(lO" 8 ) 

.584 

(10" 5 ) 

1.638 

(lO” 8 ) 

« 066 

(lO" 8 ) 

.605 

do” 5 ) 

1.623 

(lO” 8 ) 

'.067 

(lO' 8 ) 

.687 

do" 5 ) 

1.929 

H* 

O 

1 

00 

.130 

O 

1 

00 

.786 

(10“ 5 ) 

1.796 

(10" 8 ) 

.074 

(lO" 8 ) 

.909 

do" 5 ) 

1.276 

do” 8 ) 

.049 

do" 8 ) 

1.063 

do" 5 ) 

2.280 

M 

O 

1 

CD 

.349 

(10‘ 8 ) 

1.250 

(lO" 5 ) 

0.152 

do" 8 ) 

.009 

• 

(10" 8 ) 

m 


Motet $ u above has units of (ft/sec ) 2 
one foot/sec ■ 0.3C48 m/sec 




Reference meridional plane 



Fig. 2. Duct Configuration used in the Theoretical Analysis and Coordinate System 




-4 


0 


4 


8 


Fig. 3. 


kr t ' k mu r t 


Real Part of Radiation Impedence for the 
Flanged Open-end of an Annular Duct. 
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(7.62 cm) 


(10.67 cm) 


(13.72 cm) (15.24 cm) 


Characteristic Functions E (r) for an Annular 
Duot (hub-tip ratio - o.5). 
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Azimuthal angle V 0 

Sound Pressure Levels at BPF at 7 ft. (2.13 m 
from isolated Rotor due to Inflow Turbulence 





